Coherence of polaronic transport in layered metals 



o 
o 

(N 
>v 



in 

(N 



Urban LundirQ and Ross H. McKenzie 
Department of Physics, University of Queensland, Brisbane Qld 4 072, Australia 

(Dated: February 2, 2008) 

Layered systems shows anisotropic transport properties. The interlayer conductivity show a gen- 
eral temperature dependence for a wide class of materials. This can be understood if conduction 
occurs in two different channels activated at different temperatures. We show that the characteris- 
tic temperature dependence can be explained using a polaron model for the transport. The results 
show an intuitive interpretation in terms of coherent and incoherent quasi-particles within the layers. 
Further, we extract results for the magnetoresistance, thermopower, spectral function and optical 
conductivity for the model and discuss application to experiments. 

PACS numbers: 71.38.Ht,71.38.-k,72.90. 
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I. INTRODUCTION 

Layered materials show a range of interesting behav- 
ior, ranging from high temperature superconductivity to 
giant and colossal magneto-resistance. A common fea- 
ture of some of these materials (see for example Ref. 
fll l2l.{! ) is that they show a peak in the interlayer re- 
sistivity as a function of temperature. In some cases 
there is also a peak in the intralayer resistivity. We can 
identify different temperature scales, from experiment. 
T^Metermines the maxima in the interlayer resistiv- 
ity. Tj| nax determines the maxima in the intralayer resis- 
tivity. Besides this, recent angle resolved photoemission 
(ARPES) experiments^ concluded that the peak in the 
interlayer resistivity is closely related to intralayer coher- 
ence, and that there is a crossover for the spectral func- 
tion from being coherent to incoherent at a temperature 
T coh . This idea that the scattering within the layers af- 
fect the transport between the layers have recently been 
investigated by Vozmediano et al£. All of this can be ex- 



TABLE I: Temperature scales in experiment for different ma- 
terials, n.s. means that the peak is not seen in the experiment, 
indicating that, if it is there, it is higher than the temperature 
range scanned in the experiment, n.a. means the result is not 
available in our knowledge 



Material 


rr(K) 


T|P X (K) 


T coh (K) 


(Bio.sPbo.s^BaaCuaO/ 


200 


n.s. 


~180 


NaCo 2 4 a 


180 


n.s. 


~150 


Lai^Sri.eM^Oy 6 


100 


270 


n.a. 


Sr 2 Ru0 4 c 


130 


n.s. 


n.a. 


TmBa 2 Cu 3 06.4i d 


127 


n.s. 


n.a. 


(TMTSF) 2 PF 6 C 


90 


n.s. 


n.a. 


«:-(BEDT-TTF) 2 Cu(SCN) 2 / 


95 


100 


n.a. 



"from Ref. 
''from Ref. 
c from Ref. 
^from Ref. 
c from Ref. 
•'from Ref. |£ 



plained if there are two mechanisms of transport. 6,9 One, 
a coherent, dominates at low temperatures, while at ele- 



vated temperatures an incoherent contribution starts to 
dominate. In this paper we will demonstrate that po- 
laronic transport can be the mechanism providing this 
physics. This gives an intuitive explanation for the dif- 
ferent temperature scales associated with transport and 
coherence. We extend the idea presented by Alexandrov 
and BratkovskylS., and apply it to layered transport. In 
that paper they discussed (bi)polaron formation within 
giant magnetoresistance materials. 

Another powerful tool when studying polaronic trans- 
port is to study the thermopower. At high temperatures 
the conductivity is activated and the resistivity shows 
an exponential temperature dependence with a gap E a , 
the thermopower usually shows a 1/T-behavior where 
the barrier is E s . One signature of polaronic transport 
is that E s << E a , whereas for normal semiconductors 
(where the transport is activated) we have E s = Ejl^. 
By comparing the high temperature electrical resistiv- 
ity and thermopower a number of experimentalists have 
argued for the existence of small polarons in LaMn03- 
compoundsi2ii&Ii. Further evidence for small polarons 
can be found in neutron scattering data, where the po- 
laron induces a local deformation of the latticei^i 6 ^^. 
Measurements of thermopower in different directions in 
an organic quasi-two-dimensional crystal found different 
behaviour between the interlayer thermopower and the 
interlayer onei£. Further the presence of polarons was 
confirmed by photoemission experiments^. 

The approach we present is based on known approx- 
imations for the polaronsSS^i. Recent dynamical mean 
field (DMFT) calculations made on the transport of small 
polarons 22 indicate that the approximations we are go- 
ing to use overestimates the resistivity, and the exact 
functional behavior of the resistivity. The results of 
Fratini et al. 22 shows that there are two temperature 
regions. One, semiconducting region where transport is 
heavily influenced by phonon fluctuations. Then a nona- 
diabatic regime, which compares mostly to the small po- 
laron regime in the Holstein model. Here, however, we 
are more concerned with transitions between different re- 
gions of small polaron transport, not so much with the 
exact details, and it seems that the approximations we 
use does capture the essential physics. We do not claim 
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that polarons are responsible for all the observed effects, 
simply that it can provide some insight into the physics 
of layered systems. For instance, for the manganites it 
seems that the double exchange model is the preferred 
one (see Ref. |23| and references therein), although an- 
other explanation, in terms of a carrier density collapse 
due to bipolarons and their magnetic features 24 is gain- 
ing interest. Even the thermopower seems to be con- 
sistent with this model 2 ^. There are also theories using 
a combination of double exchange models and polarons 
for the localized structure 2 ^. A shorter presentation of 
some of the result from this paper has been previously 
published^. We have also investigated the problem of 
angular magnetoresistance oscillations in layered metals 
arising from incoherence 2 ^. 

The layout of the paper is as follows, in section ITT1 we 
present the model and the small polarons are introduced 
in section IIIII with the decay and Green function. In 
section IIVI we turn to the transport properties for the 
intralayer and interlayer currents and thermopower. Last 
we conclude with a calculation of the optical properties, 
section lTV Dl and a special case of the magnetoresistance 
in section HVEl 



in each layer, the only interaction between bosons in dif- 
ferent layers come from the electron tunneling, but this 
event occurs more seldom than the electron scattering of 
bosons in each layer. The Hamiltonian can be specified 
for this system. Two layers are coupled with a hopping 
Hamiltonian. Within each layer the electrons can hop 
but there is a coupling to a bosonic degree of freedom in 
each layer. We then use the Hamiltonian: 




FIG. 1: We model the two coupled layers as a anisotropic 3D 
system. Within each layer the electrons couple to bosons to 
potentially form small polarons. This is described by the the 
Hamiltonians Hi and H.2- The two layers are the coupled by 
a direct hopping term, TL t . 



II. MODEL HAMILTONIAN 



TL = TL\ + H.2 + Ti-t 



We start with a Holstein model— for an infinite system 
where the electrons interact with bosons. The Hamilto- 
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iq-Ri 



(a q + a T ), 



(1) 



'.q 



where e° is the on-site energy, w q is the dispersion of the 
bosons, tir/ is the hopping integral between neighboring 
sites i and ?y, Af q is the coupling between the bosons and 
the electrons. We want to emphasize that we will talk 
about bosons, since the theory will look the same for all 
types of bosons with a coupling given in the Hamiltonian 
above. The bosons can be phonons, spin-waves, plas- 
mons, or any other type fulfilling bosonic commutation 
rules. Since we want to study layered systems we split 
the hopping into parallel and perpendicular to the layers, 
ty and t± respectively, where i|| ^> t±. We only include 
hopping between nearest neighbors, both for the intra- 
and interlayer hopping. This enables us to write the 
Hamiltonian in a way more adapted for the layered case, 
shown in Fig. ^ The nature of the transport depends 
on how tu and tj_ compares with T, the scattering rate 
due to the bosons. We assume that that T > t±, so that 
the interlayer transport can be described by considering 
two decoupled layers. Meaning that the electrons scat- 
ters many times within the layers before jumping to the 
next2itii2L2&. We assume that we can decouple the bosons 
within each layer separately, i.e., the bosons are localized 



where 

Hi = e ° c Ui + fiu, q a q a q + *|| 4 C * + M q4 c * eiq ' Rl («< 

i q <i<q> i,q 

n 2 = J2 e °4 d i +J2 Tkj p a t a v +t \\ 4d j +Y M p4 d o ei ' l ' ni{ 

3 P <i<5> ••!> 

Tit - t±^(cjd,+h.c). 

i 

Here, and below, q, c, i,a, 1 refers to one layer, and 
p, d, j, b, 2 to the other one. 

III. SMALL POLARONS 

First we focus on the properties of the two layers sep- 
arately, i.e., we ignore the hopping term, Ti t between the 
layers. We perform a Lang-Firsov transformation^ to di- 
agonalize the Hamiltonian, excluding the hopping term, 
defined above, in each layer. Then a — > 5j = Cj-Xj, and 
dj — > dj — djYj where 



Xi = exp 
Yj = exp 



A I, 



M, 



L P 



(2) 



are the polaron operators^, for the first and second layer 
respectively. Further, a, — > dj — ^cjQ. The Hamil- 
tonian is transformed to "H = e s TLe~ s where S — 
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TSJ^ Si c l c i( a l ~ a i)- The Hamiltonian becomes 



n 



fiui, 



E^ P &j>p-E Ad H-E Ac i c 



1,7} 



3,S 



ho 



where 



A = £ 



M 2 
q _ £ o 



(4) 



is the polaron binding energy. The intralayer hopping 
terms can be treated by adding and subtracting to the 
Hamiltonian a term^I 



H S p = t\\ E ( X i X j)4 c j = E £kC k Ck 



(5) 



<ij> 



where (..) denotes a thermal average over boson states 
and this term describes a tight-binding band of small 
polarons for a square lattice within each laye r 30 i 31 



- e ' N ' lE "(^) (1+2 " B) i||[cos(fc x a) +cob(A v o)] > 

where a is the lattice constant within the layers, N 
is the number of sites in one layer, and ub(T) = 
(exp(hujq/kBT) — l) -1 is the Bose function. We see that 
the quasi-particles are described by an tight binding en- 
ergy, where the bandwidth is reduced due to the polaron 
formation. Polaron transport narrowing has been seen 
experimentally in muon-experiments^ 

There is then a residual interactional between the po- 
larons and the bosons which is described by 



H 



p—b 



= *H ^[XiXj-iXiX^Cj, 



(7) 



<i]> 



and leads to scattering of the small polarons. 



A. Decay of the quasi-particles 

Later we will need the decay, T, so we start by cal- 
culating it. We will calculate the first contribution to 
the self energy in one layer by a method similar to the 
one used by Alexandrov and Mott^i. The first non-zero 
contribution to the imaginary part of the self-energy, £, 
comes when the polaron emits one boson and absorb one 
boson. This process is shown in Fig. [21 and is induced 
by the polaron-boson scattering from Eq. (JJJ. By using 
Fermi's Golden rule we get an expression for the decay: 

r = 27T ^ I (w q - l,n q /+i;k + q - q'|7Y p _h|?i q , rt q /; k)| 2 
q,q' 

x<5(e k -e k +q-q')- (8) 




k+q-q' 



FIG. 2: Diagram describing the first contribution to the po- 
laron decay, F. The polaron emits and absorbs one boson, 
changing its momenta. 



Here, Ti p -b is the polaron-boson interaction in Eq. J7J. 
Using this Hamiltonian we get that 



{ric, 

X5 q _q'. 



1, n q /+i; k + q q'|W p _ pft |n q , n q /;k) 



q' 



q'l4 + q-q'Ck|k> 



(9) 



To simplify this to get a energy independent expression 
we use a energy independent density of states and assume 
a k-independent coupling between the electrons and the 
bosons. We only consider a single frequency loq for rea- 
sons of simplicity; it allows us to express some of our 
results in an analytical form. Then we can define the 
dimensionless coupling 
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M 

hu>o 



that will enter our equations later. We require that g > 1 
in order for small polaronic effects to be important. 43 
Using this, we get that 



Wg 2 n B {l + n B ), 



(10) 



where W — itn is the polaron bandwidth, which is sub- 
ject to narrowing due to the renormalization of the hop- 
ping t L = t«e-^ 1+2nB \ 



B. Green function in the layer 

Let us start by calculating the electron Green function 
(GF) within one layer, ignoring the coupling between the 
layers (tj_ = 0). This gives us valuable information on 
coherence of the quasi-particles, and can be compared 
to angle resolved photo-emission spectra (ARPES). Af- 
ter performing the Lang-Firsov transformation the small 
polaron GF is 

G°(k,r) = -iQ{ T )LY,j^-^HT T c^T)c\,{Q)) 



-iQ(r)e 



(e-k — iT)ir/h 
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To get the electron GF we have to convolute this GF with 
the average over two polaron operators (TXt(i)Xj/ (0)) = 
<&n>(i). This average can be decoupled and written as an 
exponential j2£LSI 

x exp i cos[q • (R, - R</)][(1 + n B )e~ iut + n B e iuJt ] 
I q 



resulting from this GF at the Fermi wave-vector as a 
function of energy. The sum over k is done by integrat- 
ing over a flat density of states. This is what is measured 
in ARPES experiments like the one in reference for 
(Bio.5Pb .5)2Ba 3 Co20 y and NaCo204. Recently, similar 
features was seen in S^RuGji^i The coherent contribu- 
tion display a peak at the Fermi energy and k-vector 
at low temperature, indicating a coherent quasi-particle. 
The peak disappears as the temperature is increased. 



After Fourier transforming the average of the polaron 
operators, giving a sum of delta functions, we will have 
a convolution 

G(k,iu n ) = -^ *(R mi ^-^)G (k>„,)e' (k - k 



/,R m ,k' 



After some algebra we come to the following expression 



R m .k' 



E 

l — — oo 



/ /[ 2 5 , E q C0S (q- R m)v /n s(l + n B )]' 
iu) n — £k' + IhujQ + iT 



Here Ii indicates a modified Bessel function of order 
Performing the summation over R m , care has to be taken 
when considering the I — term, we get the final result 
for the GF 



G(k,w) = e' 9(1+2nB) 

I hg-y/n B (l + n B ) 



1 



E 



to — £k + iT 
1 



u - e k ' 



iT 



E ■ 



Zg\Jn B {l + n B ) 



-lhuiop/2 



lu — £k' + IfkUQ + iT 



(13) 



Note that we have written the GF as a sum of a coher- 
ent and an incoherent part. This can be compared to 
the zero temperature result by Alexandrov and Mot^i. 
At T = there arc no bosons to absorb and only I > 
contributes to the GF. Also, we can compare this to the 
nonzero temperature GF by Ciuchi et ali^. The first 
line is dependent on k, thereby describing a coherent 
part. There will be a well-defined quasiparticle peak at 
lu = £k, with a spectral weight of e~ 9 ^ 1+2nB \ The second 
and third lines contains a sum over intralayer momentum 
and are therefore independent of k, they are incoherent. 
The two contributions have different temperature depen- 
dence, the coherent dominates at low temperature, and 
the incoherent at high temperature. This means that 
there is a crossover from coherent intralayer motion at 
low temperature to incoherent intralayer motion at high 
temperatures. In Fig. [5] we show the spectral function 
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FIG. 3: The quasi-particle spectral function, 

nf(e)lm[G(k.F , e)] for a electron-phonon coupling, g = 1, for 
different temperatures. The sum over k is done by integrating 
over the density of states, which we assume is flat with a 
bandwidth W = Hhujo. There are two contributions to the 
spectra function, one coherent dominating at low tempera- 
tures, and one incoherent dominating at high temperatures. 
Similar behavior have been seen experimentally-. 

crossover temperature when the contribution from the 
incoherent starts to dominate over the incoherent one. 
This will take place when 



k B T* 



2*7 ' 



(14) 



Later we will see that this relates to the temperature de- 
pendence of the interlayer conductivity in a special way. 



IV. TRANSPORT PROPERTIES 

Let us now turn to the transport properties of the lay- 
ered material which is described by the layered Hamilto- 
nian defined above. At an applied voltage V, the current 
is given by the current-current correlation function from 
the Kubo formula 20 



I„u(eV) = ^Im. 



dte" vt {TUt)jl(0)) 



(15) 
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where j is the current operator. \i and v are directions 
in the crystal. For our system, including polarons, the 
current operator for nearest neighbor hopping is 



UaJbJc) 



i,r} 



3,8 

(16) 

The first term corresponds to hopping in layer 1, the sec- 
ond term to hopping in layer 2, and the third to hopping 
between adjacent positions in the two layers. Since we 
also are going to study thermopower we give the expres- 
sion for the energy current in the same model 



The decay is calculated as usual from the imaginary part 
of the self energy, and is given in Eq. Q1U|). We use the 
tight binding approximation, Eq. ©, to get the quasi- 
particle velocity, v(k) = in one direction and get 

the conductivity, 



/ ■ e ■ e -e\ 



ffij+t) ) d ]+V d 3 Y j+r) Y 3 



t±e 



3,6 



and e is the energy of the quasi-particle. Let us separate 
the current within the layers and perpendicular to the 
layers, since they usually show a different behaviour in 
experiment. 



A. Current within the layers 

In this section we will calculate the current within the 
layers. We split the calculation into two regimes be- 
cause if we use Eq. directly in Eq. (|15|l , the result is 
too complicated to decouple, so we split the calculation 
into low and high temperatures, where different parts of 
the Hamiltonian dominates, and we can use perturbation 
theory. 



1. Low Temperatures 

At low temperatures the transport within the layers is 
coherent. If the layers are metallic we can treat them 
in a Fermi-liquid manner and use that the conductivity 
depends on the scattering rate via the scattering time, 

r = h/r, 



11 = 5? / r(k)d 2 fc. 



(18) 



ft- 



nh g 2 ns(l + n B ) 

r2ir 

I dxdy 



sin 2 (x) 



2,- 1 + cosh[/3(e + in cos(x) + tu cos(y) — /i)] 

(19) 

2. High Temperatures 



At high temperatures the polarons are localized, the 
bandwidth disappears, and the hopping, fn, is the per- 
turbation. Utilizing Eq. I|15l) for the current, we decou- 
ple the electron operators to polaron GFs in each layer, 
G = (w-A + £ + iT)- 1 . Note that there is no It- 
dependence for the polaron GFs since they are localized 
at an energy A = e — .9^0 < 0. The calculation of 
the GF in perturbation theory is described in A ppe ndix 
lAl The 4 A-operators are decoupled as in Ref. |2(J into 
diagonal (no change of boson state) and non-diagonal 
transitions (when the boson-state changes in the hop). 
The result for the non-diagonal transitions is 



(17) 



(T r At(r)A(r)At(0)A(0)) Q 



e -2g(l+2n B ) y~" 



?\Ais(l + n B ) 



For the diagonal part the four A-operators decouple and 
cancels the —1 term above when added together. Com- 
bining the two correlators and taking the imaginary part 
we convolute the two Fourier transforms similarly to 
what was done for the GF above, so that we get for the 
current: 

/ ( W ) = ^ t 2 rf 2 e -2 9 (i+2„ B ) p de 
h 11 J_ oc 2ir 

OO 

x II [W«b(1 + "b)] e - lh ^l 2 A{e + u + lhu ) 

l—~oo 

x [iif(e) — np(e + u + Ifkuo)] ■ (21) 



The conductance is obtained as usual as mi 



w=0 

The conductivity can now be plotted. The metallic, 
low temperature, part decreases with increasing temper- 
ature and the insulating, high temperature, phase takes 
over as temperature is increased. There is a peak in the 
resistivity and a crossover from coherent to incoherent 
transport, described by Eq. I|19|) an d Eq. (|21|) respec- 
tively. We did similar plots for a range of coupling con- 
stants, g, and saw that the intralayer crossover occurs at 
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a temperature given by 

fcsljn ax ~2^. (22) 

We have used the same decay, T, for both the low and 
high-temperature limits. This approximation assumes 
that the dominant part of the scattering of the carriers 
in both limits is the electron-boson coupling. The results 
are shown below in the figures below. 

B. Current perpendicular to the layers 

Let us now turn to the current perpendicular to the lay- 
ers. The current operator for an applied field in the per- 
pendicular direction in a nearest neighbor hopping model 
(from layer 1 to 2) is given in Eq. (|16J) . We assume that 
the hopping between the layers only take place between 
nearest neighbors, see Fig. ^ Then, (Ry+a — Rj) is the 
distance between the two layers, d, since (5 = 1 for near- 
est neighbor hopping. The Kubo formula, Eq. I|15f) . gives 
that we have, to second order in t±, 

I± (eV) = ^tld 2 Y I d ^ eVT ( T tc) (r)d n (r)dj (0) Cj 
x(T T Yl(r)Y j (r)Xj(Q)X jl (0))- 



We decouple the operators in the first and second layer 
respectively. This means that the Fourier transformed 
averages of the electron operators gives rise to polaron 
Green functions 



(Tct(i) Cil (0))^G?(k,ip»), 
(Td h (t)d](0))^G° 2 (p,i Pn -ioj). 



The average of the polaron-operators (X, Y) can be de- 
coupled for the two layers separately, and written as an 
exponential $(i)j22iSI as done for the GF above. Since 
the coupling in all layers are the same, M\ = M 2 , we 
can either combine the two averages of the polaron oper- 
ators into one exponential ($(t) * $(t) = (&(t)) 2 ) or keep 

-Q^them as two separate. Then we can perform the Fourier 
transform r — * u>, and if we assume that the GFs has 
an imaginary part, as above, we get the following for the 

(2fl)terlayer tunneling current 



Ix(eV) = Hf^ e -2 ff (l+2n B ) jy 00 & £ ^0 (kj e)A {K £ + eV) [f{e) _ /(£ + eV)] 



+ 



- 1 



(l 4gy/n B (l + n B ) 

OO 

Y T i +n B ) 



de 



oo k 



Y A°(k, e) ]T A° 2 (P, e + eV) [/(e) - /(e + eV)] 



<k_ 

2tt 



Y A i(K <0 E A °(P> e + eV + ltkJ ) [/(e) - /(e + eV + lhuj a 



(24) 



k belongs to the first layer, and p to the second. A\ and 
A 2 are the spectral functions for the electron GFs in each 
layer respectively, 



The index I is a combined index for the number of bosons 
emitted or absorbed in layer 1 and 2 combined 



A°(k,e) 
A° 2 (p,e) 



(e-e k ) 2 +T2' 

r 

(e-e p ) 2 +F2- 



(25) 

To illustrate this, consider what would be obtained if 
we had not combined the two exponentials, the result 
(26) would be: 
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7x(eV) oc 2tie- E -(^) 2(1+2 " fl) - E p(^) 2(1+2nB) 



n e 

q.p (,i'=— oo 



Z 100 cfc 

; / — A 2 (e - Z'/iu; q )A l (e + eV" + Zfrw p ) [n F (e - f/iWq) - n F (e + eV + lhio p )} 
J-oo 2tt 



(27) 



Z belongs to the first layer and counts the number of 
bosons attached to the electron, V refers in a similar fash- 
ion to the second layer. The connection between Eq. (|24H . 
and Eq. (|27|l can be found from the identity for the Bessel 
functions 



From Eq. (|24|l we can extract the conductivity by a 



00 



l,l'=-c 



Ii-v{x 1 )I v (x 2 ) 



Ii(x 1 +x 2 ). (28) 



l=-c 



The expression for the current, Eq. I|24l) . has a con- 
tribution from coherent and two from incoherent trans- 
port. Note the similarity in structure of Eq. (|24(l to the 
expression for the GF, Eq. The first term corre- 

sponds to transport which conserves the intralayer mo- 
mentum in the tunneling process. This is seen since the 
crystal momentum k is the same for the spectral func- 
tion for the two different layers. For the other terms, 
the intralayer momentum is not conserved, (in each layer 
the sums over the momentum are separate). The sec- 
ond row corresponds to transport when the net number 
of bosons in the system is unchanged. When the quasi- 
particle tunnels it leaves behind the cloud of bosons in 
one layer and attaches to a replica of bosons in the sec- 
ond layer. The third row describes transport when a net 
number of bosons is absorbed (I > 0) or emitted (/ < 0), 
thus changing the energy of the polaron in the hop be- 
tween the two layers. In a recent paper™ we established 
a connection between the intralayer coherence and the 
appearance of dips in the angular magnetoresistance, the 
so called "magic angles". We showed that a contribution 
from incoherent jumps between highly conducting one- 
dimensional strands of molecules gives a natural expla- 
nation of the phenomena observed in the magnetoresis- 
tance. At low temperature the coherent part dominates 
but at high temperature (high compared to the boson 
energy, Hljq, ksT > Hujq) the incoherent mechanism of 
transport will dominate. Thus, there is a crossover from 
coherent to incoherent transport. The crossover temper- 
ature is fixed by having equal contribution from the co- 
herent and the incoherent contributions. Ignoring the 
contribution from the I ^ terms in Eq. i|24|) we can get 
an approximate expression for the crossover temperature 
as: 



2 3 5 



1. 68 ^2. (29) 



simple derivative mi = e ,f 7 JU . In Fig. 0] we plot 

11 a\ev ) eV=Q 

the conductivity as a function of temperature for one 
value of g. The crossover is clearly seen. 




FIG. 4: Interlayer resistivity as a function of temperature for 
different values of the coupling g . At low temperatures, the 
transport is predominantly coherent, as seen from Eq. iTHl . 
Then, as the temperature is increased, the incoherent mecha- 
nism of transport take over. We have crj_ = a c f h ' + u™ c ', the 
two contributions are shown separately and together in the 
plot. The crossover from coherent to incoherent transport is 
clearly seen, 1 = ^-t^d 2 , W = 77hLu . 



In general, the interlayer conductivity for identical de- 
coupled layers isS, 



2e 2 , 

a - = X* 



J de£|A(k,e) 



de 



(30) 



where A(k, e) is the electron spectral function for a single 
layer. Directly substituting Eq. (|13H in Eq. (|30fl we obtain 
the same result found from Eq. (|24|) . 
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We can check the result by taking some limits in 
Eq. H24|) When the temperature is zero, we get 

a ±i T = o) = 



2e^ 



-il e 



2r 



(31) 



Thus, at low temperature when only the first (coherent) 
term contributes to the conductivity the temperature de- 
pendence of <7j_ is governed by the temperature depen- 
dence of the decay, T, given by Eq. (jTU|l . for the polaron 
case. If we expand Eq. I|24|) for high temperatures the 
conductivity behaves approximately as: 



(7j_ oc 



y-3/2 



(32) 



which is consistent with the equipartition theorem^. 
If we take the limit g = we get: 



9 ^ Z" 00 ri 

a ± (g = 0) = ^-iid 2 £ / ^(k^fre) 

x/3n F (e)[l-n F (e)], (33) 
as expected from transport theory^. 

C. Thermopower 

Let us now turn to calculating the thermopower for 
intralayer and interlayer transport. The thermopower is 
defined as a correlator, using the heat current instead of 
the electrical current in Eq. (|15|) . (See Ref. l2(if) 

and, using that the current-current correlator gives us 
L n , and a = L n /T, we get 



S 



T L 11 ~ T 2 a 



(35) 



For the intralayer thermopower we consider the low and 
high temperature limits separately. 



1. Low temperature intralayer thermopower 

At low temperatures the correlator is similar to the one 
calculated for the intralayer low-temperature conductiv- 
ity, except for an additional efc/2 in the (energy)current 
operator. This factor only contributes an additional 
iy cos(fc;ca)/2 if we assume that we do the measurement 
along x. The result is that, for low temperatures, the 
thermopower is 

clow 1 H I d 2 kf{ekk) [1 - f{e kk )]sm 2 (k x a)cos(k x a) 



Tie 



J d 2 kf(e kk ) [1 - f(e k k)] sm 2 (k x a) 



(36) 



2. High temperature intralayer thermopower 



At high temperatures, we can follow the same steps as 
for the intralayer current with the only difference that the 
energy operator is the current operator multiplied by ^jf- . 
Then, in Eq. the correlators cancel, and we simply 
end up witbS: 



2eT' 



(37) 



The two results in the low and high temperature re- 
gions, Eq. and Eq. (|37|l respectively, both falls off 
as t^. This means that there would be no peak in the 
intralayer thermopower corresponding to any transition 
between coherent and incoherent transport. Therefore, 
the transition in the intralayer transport is more clearly 
seen in the electrical transport, not the thermopower. 
The 1/T-dependence is typical for polarons at high tem- 
peratures as seen, e.g., in La 2 /3Cai/ 3 Mn03 flimsy and 
(La,Ca)Mn0 3 — 



3. Interlayer thermopower 

We follow the same steps as for the interlayer conduc- 
tivity with the replacement of one current operator by 
one energy-current operator as in Eq. I|34() . The result 
for L 12 is 
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12 



2e 



t 2 e -2 9 (l+2n B ) d 2 



{/"OO i 
/ _ ^ E ^^?(k, e)A°(k, e + eV) [/(e) - /(e + e7)] 



(j [4. 9v /n B (l + n B )] -l) /" ^?(k, e) £ A°(p, e + eV) [/(e) - /(e + eV)] 



OO 



rfe 



— V £ kJ 4°(k, e) V 4°(p, e + e y + Z^o) [/(e) - /(e + eV + Z^ ) 

Z7T * — ' * — ' 



(38) 




FIG. 5: The top panel shows the interlayer resistivity as 
a function of temperature for different electron-phonon cou- 
pling strengths, p^ 1 = ^-t^d 2 , W = 160fkj . The peak 
correspond to the transition between coherent and incoherent 
transport. The lower panel shows the interlayer thermopower, 
S° = . The peak in the thermopower occurs at a lower tem- 
perature than for the resistivity. The inset shows the ther- 
mopower for small electron-phonon coupling, note that it can 
change sign although we do not consider carriers of hole type 
here. At high temperatures, and strong electron-phonon cou- 
pling, the thermopower decays exponentially. 



Here £k = ek — ji. The thermopower is then given by 
Eq. 1351) . In Fig. [S] we make a comparative plot of the 
resistivity and the thermopower between the layers. For a 
Fermi-liquid the thermopower would fall off as y at high 
temperatures (see, e.g., Salamon et al>al), fitting a curve 
to our numerical results shows that our expression for the 
interlayer thermopower falls off exponentially instead. 



D. Optical conductivity 

The optical conductivity is given by calculating the 
derivative of the frequency dependent current in 
Eq. I|24l) with respect to u> = eV. Assuming that the 
following relations hold (when T <C W) 



J dxD(x)A(x,e) = J dxD{x)-^ 



(e-x) 2 +T 2 
dxD(x)A(x, e)A(x, e + lu) = 

r r 



= D(e) 



dxD(x) 



(e-x) 2 + T 2 (e + to-x) 2 +T 2 



D(e)- 



,2 +T 2 



We get that the full expression for the optical conductiv- 
ity is: 
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h 



d r 

27T + 1 J 



+ (/o [4 5 V«s(l + n B )] - l) y ^£>(e)£>(e + Lu)/3n F (e + u) [I - n F (e + u)} 
+ y h 4 9v /n B (l +n B ) e - lhuJol3/2 / — D(e)D(e + u + l?uJo)f3n F (e + u + lfrLJ )[l-n F (e + u + lHuj 

2r w 



i= — oo 
f0 ° de 



— OO 



[n,F(e) - n f (e + u)] 
de 



+ (la [^9y/n B (l + n B )\ - l) y ^ J D(e)D'(e + w) [n F (e) - n F {e + u)] 

°° r i f°° de 

+ Y h 4. 9v /n B (l+n B ) e ~^ 0/3 / 2 / — D{e)D\e + u> + ltlw ) [n F {e) - n F {e + u> + IHlj )} 

I 



Here D' is the derivative of the density of states (which 
is zero for a constant DOS). In the figure for the op- 
tical conductivity, Fig. [U we see that there is a large 
coherent Drude peak at low frequency, which disap- 
pears as the temperature increases, consistent with ex- 

• ji _ * i 3 mi_ _ j; 




0.3 
"S 0.2 

e 
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Incoherent 
~ v — 

\ 

1 - 
1 







0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 0.5 

co/4 t|| 

FIG. 6: The top panel shows the frequency dependence of 
the interlayer optical conductivity for two different tempera- 
tures. The Drude peak (at uj — 0) disappears when coherence 
is lost, due to the destruction of coherent quasi-particles with 
increasing temperature. The two lower panels shows the opti- 
cal conductivity divided into the two contributions, coherent 
and incoherent plotted for two different temperatures. The 
input density of state is flat with a bandwidth of W = 77huo- 

ance of the Drude peak is consistent with the peak 
in the spectral function disappearing, Fig. [3] Such a 
behavior has been observed^ in Sr2Ru04, where the 
Drude peak disappeared above 100K. Another system 
is Lao.825Srg.i75Mn04 where the Drude peak disappears 



above 200Kia 



E. Interlayer Magnetoresistance for a field parallel 
to the layers 

We can also make a statement about the magnetoresis- 
tance in a certain limit. If we apply a magnetic field, B, 
parallel to the layers (the x-y plane) there is an orbital ef- 
fect on the paths of the electrons. This can be described 
by a shift in wave vector, k — > k -|A, where A is the 
vector potential for the magnetic field. For a magnetic 
field in the x direction, when an electron tunnels between 
adjacent layers it undergoes a shift in the y-component 
of its wave vector by —dB 3 %<. In the general expression 
Eq. I^n)) \A(k, e)| 2 is replaced with an equation contain- 
ing Ai(\c, e)^2(k+ j^dBy), since there will be a difference 
in the vector potential between the two layers. 

However, since the incoherent part of the conductivity 
contains a summation over k-space and is independent 
of k, this will be unaffected by the magnetic field, i.e. 
E P 4KP + JidBy, e + eV) = £ p A°(p, e + eV) since the 
sum span over the first Brillouin zone. 

Thus, we will have two contributions to the interlayer 
conductivity and one is ^-independent: 

<t_l(B) = a c f h {B) + < coh (B = 0). (39) 

a coh (B) decreases with increasing magnetic fieldS^i 



y/l + {evpcBTf 



(40) 



where v F is the Fermi velocity. If we increase B, the co- 
herent part decreases, and, therefore, T™ ax would shift 
to lower values. A separation of the conductivity in two 
parts, as in Eq. J2HJ), has been proposed previously^ on 
a phenomenological basis, in order to describe the mag- 
netoresistance of Sr2Ru04 (Except there a weak field de- 
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TABLE II: Temperatures scales obtained in the small polaron 
model applied to layered systems. 

fcoh ~ o"incoh GF coherence Interlayer transport Intralayer transport 



pendence is associated with the incoherent contribution 
due to Zeeman splitting). 

V. DISCUSSION 

We have presented a layered polaron model for sys- 
tems consisting of two-dimensional layers coupled by tun- 
neling. We have found that when the temperature is 
lower than the characteristic boson frequency the physics 
is dominated by coherent transport where the electrons 
scatters of bosons in the layers. Upon increasing the 
temperature a transition is made into a region where the 
physics is governed by incoherent small polarons. The 
small polarons are localized at the lattice sites and hop 
to new sites. We have extracted results for intralayer 
and interlayer transport, thermopower, ARPES, optical 
conductivity and magnetoresistance. 

In Fig.|7|we plot the different crossover temperatures as 
a function of the dimensionless electron-boson coupling. 

T^tin rlifl-nvnTnf tpmnprnfuro or>cil^c Qccn^idtn/l inrifti email 



FIG. 7: Crossover temperatures as a function of electron- 
phonon coupling constant, for the intralayer and interlayer 
resistivity, and the interlayer thermopower, For a fixed cou- 
pling and increasing temperature first the coherence between 
the layers is lost, then there is a crossover peak in the inter- 
layer thermopower, and last the coherence within each layer 
is lost at elevated temperatures. The sequence of crossover 
does not change if t\\jW is changed. 



A- -AT 



O T Thermopower 



O- O O O O O-..0. 



❖ . 



O O.. 



° ° o o o O O < 



O- O- O - o- 



polaron transport is summarized in the table ITTl 



The presented theory differs in the way polarons are 
formed from the theory by Ho and Schofield 4 -. In their 
paper they assume that the polarons are formed from 
phonons traveling along the c-direction. Some of the re- 
sults are similar, but in principle it is a different theory. 
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APPENDIX A: GREEN FUNCTION IN THE 
HIGH TEMPERATURE LIMIT 

In order to calculate the intralayer conductivity at high 
temperatures, when polarons are formed and are local- 
ized, we need the GF. We perform perturbation theory 
in t and include the diagrams shown in Fig. |H1 When 




FIG. 9: Example when the polaron can hop 3 steps away, 
from the solid circle to the site marked with a box. From the 
sites marked with a there are 3 paths back (marked 1-3) and 
there are 3 paths there, there are 8 such equivalent a-sites. 
The polaron can also take the path to the site marked /3, there 
are 4 equivalent sites of this type. 



FIG. 8: The diagrams taken into account when calculating 
the self-energy for the trapped polaron in the high tempera- 
ture limit. 



possible paths are 4+2-4n(n-l) 
self energy is 



Thus we get that the 



summing the series shown in Fig. |SJ we have to consider 
the lattice the electron moves in, i.e. the number of possi- 
ble ways to return to the original site. In Fig. we have 
shown the number of possible paths when the polaron 
hops 3 jumps away from its original site. We have to find 
a general expression for the number of possible jumps for 
the whole lattice, and take into account identical paths. 
Generally if the electron hops 2n steps, the number of 

o — o — o — -o — o — o — o 



a 



P 



a 



o — o — o — o — o — o 



o — o — o — o — o- — o 



a 



a 



(3 



o — o — o- -o- — o — o 




o 
o 
o 



E(w,r) = 4^[l + 2n(l-n)]i 2 "G 2n - 1 



= At 2 G 



1 + t 2 G 2 + i 4 G 4 
(l-i 2 G 2 ) 3 ' 



(Al) 



where G is the local polaron GF, and t is the hopping 
integral within the layer. 
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